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Abstract 

Background: Mycobacterium avium subsp. avium {Maa) and M. avium subsp. hominissuis {Mah) are environmental 
mycobacteria and significant opportunistic pathogens. Mycobacterium avium infections in humans and pigs are 
mainly due to Mah. It is not known whether this is caused by a difference in virulence or difference in exposure to 
the two subspecies. The aim of the present study was to investigate the ability of the M. avium subspecies to 
replicate intracellular^ and to characterise the gene expression program triggered by infection of human primary 
macrophages. 

Results: All isolates were able to invade and persist within human macrophages. However, intracellular replication 
was only evident in cells infected with the two Maa isolates. Transcriptional responses to the isolates were characterized 
by upregulation of genes involved in apoptosis, immune- and inflammatory response, signal transduction and NF-kB 
signaling, cell proliferation and T-cell activation. Although similar pathways and networks were perturbed by the different 
isolates, the response to the Maa subspecies was exaggerated, and there was evidence of increased activation of type I 
and II interferon signaling pathways. 

Conclusion: Mycobacterium avium isolates of different genetic characteristics invaded monocytes and induced different 
degree of macrophage activation. Isolates of Maa were able to replicate intracellular^ suggesting that differences in 
exposure, uptake or induction of adaptive immunity are more likely explanations for the difference in prevalence between 
M. avium subspecies. 
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Background 

The Mycobacterium avium complex comprises a group 
of environmental and opportunistic intracellular patho- 
gens [1] with the ability to persist within macrophages 
and escaping the host's killing mechanisms [2]. The 
complex is constituted by the subspecies of M. avium 
together with M. intracellulare [3]. Two subspecies of 
clinical importance for humans and animals are M. 
avium subsp. hominissuis {Mah) and M. avium subsp. 
avium {Maa) [4]. Mah is the subspecies mainly isolated 
from porcine and human cases of MAC infection [4], 
while Maa is the causative agent of avian tuberculosis, 
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but is occasionally responsible for M. avium infections in 
mammals [5]. Little is known about the routes of trans- 
mission and the virulence mechanisms of M. avium, and a 
zoonotic potential has not been ruled out [6]. 

Definite reasons for Mah being the dominating sub- 
species in humans and pigs infected with M. avium have 
not been identified, although differences in exposure 
could be considered a contributing factor. Explanations 
to the higher prevalence of Mah in the human and 
porcine population might also be found in differences 
between Mah and Maa at the molecular level of inter- 
action between the host and mycobacterium, a process 
crucial to the outcome of infection, which therefore 
might affect the prevalence of the subspecies. However, 
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little has been done to investigate whether Mah and 
Maa interact differently with their target cells. 

Mycobacterial virulence depends on the ability of these 
bacteria to invade, persist and replicate within the hostile 
environment in macrophages of the host organism [7]. 
This is achieved by the protective properties of the my- 
cobacterial cell wall and by interference with the host's 
immune response [8,9]. To assess the ability of Mah and 
Maa to grow intracellularly, the use of primary cells will 
give a more representative result compared to immorta- 
lised macrophage cell-lines. Tissue resident macrophages 
represent a heterogeneous population and it is not clear 
whether M. avium infect a particular phenotype. Mono- 
cyte derived macrophages matured in absence of any stim- 
uli give rise to a mixed culture of both the Ml and M2 
phenotype [10] and such cultures may thus be considered 
appropriate for infection studies with M. avium. Finally, 
using macrophages derived from the relevant target spe- 
cies is important since the global inflammatory response 
evoked in mice recently was demonstrated to be vastly 
different from humans [11]. To the authors' best know- 
ledge, comparative studies of the ability of Maa and Mah 
to invade, replicate and trigger immune responses in 
human primary macrophages have not been published. 
The aim of the present study was therefore to assess the 
growth of M. avium isolates in human primary macro- 
phages and to characterise the response triggered in in- 
fected macrophages by selected isolates. The isolates used 
in the present study represented the two different subspe- 
cies, were derived from different host species and differed 
in biofilm abilities, glycopeptidolipid (GPL) genes, and 
in the presence of the insertion sequence (IS) element 
ISMpal [12]. The present data showed that Maa repli- 
cated to a larger extent in human macrophages suggesting 
that an inability to grow intracellularly is not likely the 
cause of the discrepancy of the prevalence of the two sub- 
species in humans. Furthermore, the response induced in 
the macrophages could not entirely explain the difference 
in growth rate. 

Methods 

Mycobacterial isolates 

Four clinical isolates of M. avium subsp. hominissuis 
{Mah VI101, Mah 1655, Mah HI and Mah H38) and 
two clinical isolates of M. avium subsp. avium (Maa 
1794 and Maa 1553) were used in the present study. 
Details on origin, biofilm producing abilities, presence of 
ISMpal, nsGPL- and ssGPL genes were described by 
Johansen et al. [12], and are summarised in Table 1. The 
human isolates were from a strain collection at the The 
Norwegian Institute of Public Health. All isolates were 
grown in Middlebrook 7H9 (BD Diagnostics, Sparks 
MD) w/OADC (BD Diagnostics), 1% Tween 80 and 10% 
glycerol under agitation for two weeks before aliquots 



Table 1 Characteristics of Mycobacterium avium isolates 



used in the present study 


Isolate 


Subspecies 


Origin Biofilm 


\SMpa1 


nsGPL 


ssGPL 


VI101 


Hominissuis 


Porcine + 




+ 




1655 


Hominissuis 


Porcine 


+ 






HI 


Hominissuis 


Human 




+ 




H38 


Hominissuis 


Human 




+ 


+ 


1553 


Avium 


Avian 




+ 


+ 


1794 


Avium 


Avian 




+ 


+ 



Isolates of Mycobacterium avium used in the present study show different 
combinations of subspecies designation, origin, biofilm abilities and the genomic 
presence of ISMpal and genes involved in the synthesis of non-specific (ns) or 
serovar specific (ss) GPLs [12]. 



were stored at -70°C in an equivalent medium contain- 
ing 10% glycerol. The concentration of bacteria in the 
stock suspensions was determined by colony forming 
unit (CFU) counts. Bacterial inocula for the gene expres- 
sion assay were prepared from seven days subcultures 
grown on plates of Middlebrook 7H10 (BD Diagnostics) 
with 10% OADC (BD Diagnostics) at 37°C. To adjust the 
number of bacteria in the inoculum, hsp65 real-time 
PCR [13] was performed on serial dilutions in HBSS 
(Invitrogen, Oslo, Norway) of a suspension adjusted to 
McFarland standard 2.0, homogenised through a 23 G 
needle to minimise clumping. 

Cell culture 

Peripheral blood mononuclear cells (PBMCs) were iso- 
lated from buffy coat drawn from healthy human blood 
donors. The material was commercially obtained from a 
blood donation centre, and the study was approved by 
the Regional Committee for Medical Research Ethics, 
South-East Norway. Fifty mL buffy coat was diluted 1:1 
with RPMI 1640 (Sigma-Aldrich, Oslo, Norway) with 
1% L-glutamine (Sigma-Aldrich). PBMCs were harvested 
by density gradient centrifugation with Lymphoprep 
(Medinor, Bryn, Norway) and resuspended in PBS with 
0.5% BSA and 2 mM EDTA at pH 7.2. Isolation of CD14+ 
cells was performed using CD 14 MACS magnetic beads 
(Milteny Biotec, Bergisch Gladbach, Germany) according 
to the manufacturer's instructions. The cells were resus- 
pended in RPMI 1640 with 1% L-glutamine supplemented 
with 10% human serum to the desired concentration of 
the respective infection assay, added to the appropriate 
Costar 8 cell culture plates (Sigma-Aldrich) and incubated 
at 37°C with 5% CO2. The CD14+ cells were incubated 
overnight before used in the uptake and replication assays 
and allowed to mature into macrophages by incubation 
for five days before used in the gene expression study. 

Uptake and intracellular replication 

Uptake and intracellular replication of the six M. avium 
isolates in human blood derived CD 14+ cells were assessed. 
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The assays were performed in 96-well cell culture plates 
(4 x 10 s cells/mL, 125 ul per well), which in our hands gave 
a stable monolayer for most donors. Each isolate was tested 
in triplicates at an MOI of 5:1. The cells were washed twice 
with pre-warmed medium to remove extracellular bacteria 
before the quantification of intracellular bacteria was per- 
formed. Intracellular uptake of the bacterial isolates was 
assayed after three, six, 12 and 24 hours in cells from three 
donors. 

Intracellular replication of the same isolates was assayed 
in six other donors, following the protocol described above. 
Based on the results from the uptake studies, CD 14+ cells 
were incubated with the bacterial isolates for six hours. 
After removal of extracellular bacteria, fresh medium was 
added and the cells were incubated further. Cell lysates 
were harvested at six hours, and after one, four and seven 
days of infection. 

Quantitation of bacteria and cells 

Numbers of intracellular bacteria and cells was deter- 
mined by enumeration of mycobacterial and human ge- 
nomes in cell lysates. Replication of bacterial isolates in 
the supernatant had been tested and found insignificant. 
Lysis of cells was achieved through adding milliQ water 
to the wells. After 30 minutes, lysates were transferred 
into 2 mL O-ring vials (Biospec products, Techtumlab, 
Umea, Sweden) containing 200 uL silica beads (Biospec 
products) and stored at -20°C. Lysates from triplicate 
wells were pooled. After completing all sampling, ly- 
sates were thawed and inactivated at 100°C for 30 min, 
prior to mechanical disruption using the MiniBead- 
Beater (Techtumlab) for 45 sec. 

Enumeration of mycobacteria and cells was achieved 
by performing a duplex real time PCR of a 103 bp long 
segment of the mycobacterial hsp65 gene and a 150 bp 
long segment of the human p-globin gene on the lysates, 
as described by Salte et al. [13]. Briefly, the segment of the 
mycobacterial gene hsp65 (GroEl2) was amplified using 
the primers MycoFPl (5 CGAGGCGATGGACAAGGT- 
3') and TB 12 (5 CTTGTCGAACCGCATACCCT-3 '), 
and the VIC tagged TaqMan-MGB probe MycoPrl (5'- 
AACGAGGGCGTCACCGTCG-3 ') [14]. The segment of 
the human p-globin gene was amplified by using the 
primers BG-F (5 '-TGCCTATCAGAAAGTGGTGGCT- 
3') and BG-R (5 -GCTCAAGGCCCTTCATAATATCC- 
3'), and the FAM tagged TaqMan-MGB probe BG-TAQ 
(5 -TGGCTAATGCCCTGGCCCACAA-3 ') [13,15]. Dilu- 
tions of known concentrations of commercially available 
mycobacterial DNA (ATCC-19015D-5, LGC Standards, 
Middlesex, UK) and human DNA (Applied Biosystems, 
Foster City, CA, USA) were run together with the samples 
and used to create standard curves. The conversion of the 
known amount of standard DNA to number of bacterial 
and human genomes were described by Salte et al. [13], as 



well as the validation of the number of mycobacterial tar- 
gets against CFU counts. Total volume of each reaction 
mixture was 20 ul, consisting of 8 ul lysate, 10 ul Quanti- 
tect" mastermix (Quiagen, West Sussex, UK) and 1 ul of a 
primer/probe mix for each target, with a stock concentra- 
tion of 8 uM for each primer and 4 uM for the probe, 
resulting in a final primer- and probe concentration of 
0.4 uM and 0.2 uM, respectively. Samples and standard 
DNA dilutions were run on a Stratagene Mx3005P MxPro 
v 4.10 (Stratagene, La Jolla, CA, USA) and data analysed by 
the Stratagene MxPro v 4.10 software (Stratagene) with the 
following thermal cycling conditions: 95°C for 15 min, 
followed by 40 cycles of 95°C for 1 min and 60°C for 1 min. 

Gene expression profiling 

Isolates Mah VI101, Mah 1655 and Maa 1794, were 
used to infect macrophages for gene expression analysis. 
CD 14+ cells were obtained as described above and 
allowed to mature into macrophages by incubation for 
five days in duplicate wells (6 well plates, 1 x 106 cells/ 
mL, 2 ml/well). The cells were then infected with the 
three isolates at a MOI of 10:1 for four hours. Total 
RNA was isolated of the adherent macrophages using 
the RNeasy Mini Kit 8 (Qiagen, Hilden, Germany), fol- 
lowing the manufacturer's instructions. RNA from the 
different donors within each assay was pooled per bac- 
terial isolate. The first assay was performed in macro- 
phages from three donors and subsequently repeated 
twice in two new donors per assay. The quality and 
integrity of RNA was confirmed by Agilent 2100 Bioana- 
lyzer (Agilent Technologies) and Nanodrop™ (Thermo 
Fisher Scientific, Wilmington, DE, USA). Gene expres- 
sion levels were profiled using Illumina HumanHT-12 
Expression Bead Chip (Illumina Inc., San Diego, CA, USA). 
The gene expression data was exported from Illumina's 
GenomeStudio software package without background 
subtraction or normalization. The raw data was prepro- 
cessed in the R language and environment for statistical 
computing [16] employing the normexp by control 
(neqc) algorithm [17]. Probes that were not significantly 
expressed above background (detection p-value < 0.01 
in at least one sample) were filtered out of the analysis. 
Differentially expressed genes were identified employing 
Bayesian/moderated i-statistics with False Discovery 
Rate (FDR) control for multiple testing [18]. Differen- 
tially expressed genes were identified at FDR < 0.05. The 
raw microarray data have been made publicly available 
in the Array Express repository (http://www.ebi.ac.uk/ 
arrayexpress/ experiments/E-MTAB- 1 101 /) . 

Bioinformatics 

Biological functions enriched in the differentially expressed 
genes were identified using modular enrichment analysis 
with DAVID (Database for Annotation, Visualization and 
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Integrated Discovery) [19,20]. Enrichment scores > 1.3 
were considered statistically significant. Biological path- 
ways enriched in the data were identified using IPA soft- 
ware (Ingenuity Systems, Redwood City, CA). Upstream 
regulator analysis was performed in IPA software to 
infer the putative pathways that give rise to the ob- 
served gene expression changes. Molecular interaction 
networks were reconstructed from the differentially 
expressed genes in IPA software using mechanistic data 
from prior studies [21]. In this analysis differentially 
expressed genes were defined at FDR < 0.01. 

Protein analysis 

CD 14+ cells were isolated and matured as described above 
for gene expression studies. Cell lawns were infected at an 
MOI of 5:1 for 24 hours with Maa 1794, Mah VI101 and 
Mah 1655. Supernatants were removed for cytokine pro- 
filing, and subsequently the adherent cells were washed 
twice in PBS and collected with a cell scraper in 100 ul of 
lx SDS Sample Buffer for immunoblotting. All samples 
were stored at -70°C until further processing. 



Results 

Uptake and intracellular replication 

All the isolates were taken up by the CD14+ cells, and 
intracellular bacteria were detected after 3 hours (Figure 1). 
Maa 1553, Mah H38 and Maa 1794 appeared to be taken 
up to a greater extent, while the Mah isolate VI101 showed 
the lowest degree of uptake in macrophages. Between six 
and 12 hours of infection some of the isolates showed an 
exponential increase in the number of intracellular bacteria, 
suggesting the event of onset of intracellular replication. Six 
hours was therefore chosen as the time between inocula- 
tion and removal of extracellular bacteria for the subse- 
quent study of intracellular replication over time. At this 
point in time VI101 had significantiy (p < 0.05) lower num- 
bers of intracellular bacteria compared to the other isolates 
except for Mah HI, while no significant difference were ob- 
served between the other isolates using the paired t-test. 

Intracellular growth was only observed for the two 
Maa isolates. There was a significant increase in the ra- 
tio of mycobacteria to human genomes at day four for 
isolate Maa 1794 (p = 0.04) and at day 7 for Maa 1794 
(p = 0.01) and Maa 1553 (p = 0.02) and in absolute 



SDS-PAGE and immunoblotting 

Proteins were separated by horizontal SDS-PAGE with 
precast gradient ExcelGel™ SDS Gradient 8 - 18% (GE 
Healthcare Life Sciences, Uppsala, Sweden) and ExcelGel™ 
SDS Buffer strips (GE Healthcare) in a Multiphor " II 2117 
(GE Healthcare), prior to diffusion blotting onto nitro- 
cellulose membranes with pore size 0.2 um (Schleicher & 
Schuell, Dassel Germany) [22]. Proteins remaining 
within the gel were visualized by Coomassie brilliant 
blue staining. Membranes were blocked with PBS with 
2% BSA and incubated over night with primary anti- 
bodies (Cell Signalling Technology, Danvers, MA) rec- 
ognizing the following selected proteins involved in 
apoptosis: Cleaved caspase-3 (Aspl75), cleaved caspase-8 
(Asp391), cleaved caspase-9 (Asp315), cleaved caspase-9 
(Asp330), Bcl-2, Phospho-Bcl-2 (Ser70). Detection of 
bound antibodies was performed by incubaction with 
ECL T " Anti-rabbit IgG, Horseradish peroxidise labelled 
linked F(ab')2 fragment from donkey (GE Healthcare) and 
visualization by Immun-Star™ WesternC™ Kit (Bio-Rad, 
Hercules, CA). 



Cytokine profiling 

The concentration of TNF-a, IL-23(pl9), IL-10, IL-8 and 
IL-6 was measured in supernatants by using Bio-Plex™ 
Cytokine Assay (Bio-Rad) within the Bio-Plex suspension 
array system (Bio-Rad), following the manufacturer's in- 
structions. Differences between isolates were tested for sig- 
nificance using the Wilcoxon matched pairs signed-rank 
test, considering p-values < 0.05 significant. 
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Figure 1 Uptake of six clinical isolates of Mycobacterium avium 
in human CD14+ cells. Cells were inoculated for three, six, 12 and 
24 hours. Numbers represent the mean values from three donors, 
showing the ratio of numbers of mycobacterial to human genomes 
(A) and absolute numbers of mycobacteria per well (B), measured 
by real-time PCR of single copy genes in cell lysates. 
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numbers of bacteria for isolate Maa 1794 (p = 0.01) at 
day 7 using the paired t-test (Figure 2). The remaining 
isolates, which were all Mah, did not show significant in- 
crease in intracellular numbers throughout the duration 
of infection, however the count of intracellular bacteria 
was quite stable, indicating that the bacteria were able to 
persist within the macrophages. 

Enumeration of both cells and intracellular mycobac- 
teria was performed by real-time PCR of single-copy 
genes. This method does not distinguish between viable 
and non-viable bacteria and cells, but it has the advan- 
tage that the number of bacteria is not influenced by the 
degree of clumping which varies between Maa and 
Mah. Furthermore it enabled the measurement of the 
number of intracellular bacteria per cell in addition to 
absolute counts of bacteria. Some degree of cell death 
and loss of adherence were observed in both infected 
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Figure 2 Intracellular replication of six clinical isolates of 
Mycobacterium avium in human CD14+ cells. Sampling was 
performed after six hours of inoculation with added bacteria, then 
one, four and seven days after removal of extracellular bacteria. 
Numbers represent the mean values from six donors, showing the 
ratio of numbers of mycobacterial to human genomes (A) and 
absolute numbers of mycobacteria per well (B), measured by real-time 
PCR of single copy genes in cell lysates. There was a significant increase 
in the ratio of mycobacteria at day four for isolate Maa 1 794 (p = 0.04) 
and at day 7 for Maa1794 (p = 0.01) and Maa 1553 (p = 0.02) and in 
absolute numbers of bacteria for isolate Maa 1 794 (p = 0.01) at day 7 
using the paired t-test. 



and non-infected cells, and varied considerably between 
donors. Similar results were observed whether absolute 
numbers of bacteria or the bacteria per cell ratio were 
used, however the difference in cell death between do- 
nors may explain why the absolute counts of bacteria 
only reached a significant level of p < 0.05 for Maa 1794 
at day 7, while the ratio was significant for both Maa 
isolates at day 7 and for Maa 1794 at day 4. 

Gene expression analysis and bioinformatics 

Based on the results from the initial replication assays, 
three isolates were selected for gene expression studies 
in macrophages. Maa 1794, which replicated intracel- 
lularly, and Mah VI101, which only persisted within 
cells, were selected, in addition to Mah 1655, which, as 
opposed to the other two isolates, contains the IS elem- 
ent ISMpal and belong to an RFLP type that has been 
demonstrated to be particularly virulent in pigs [23,24]. 

The Illumina 8 HumanHT-12 v3 Expression Bead al- 
lows whole-genome expression analysis of 12 samples 
on the same chip, each array targeting genomewide with 
more than 48,000 probes. Of the 14,581 probe sets with 
signal intensity above background, 2,766 probe sets were 
differentially expressed (FDR < 0.05) in the cells infected 
with M. avium (1,229 up-regulated, 1,537 down-regulated), 
regardless of the individual isolate, compared to uninfected 
control cells (Additional file 1: Table SI). The largest 
response was induced by Maa 1794 where 2,439 probe 
sets (representing 2,147 unique genes) were differentially 
expressed. Of these probe sets, 1,077 were up-regulated, 
and 1,362 were down-regulated. Mah VI101 modulated 
expression of 1,520 probe sets (810 up-regulated, 710 
down-regulated) representing 1,348 unique genes. In com- 
parison, Mah 1665 altered the expression of only 341 
probe sets (239 up-regulated 102 down-regulated) repre- 
senting 298 unique genes. A large overlap was observed in 
the responses elicited by the three isolates (Figure 3). There 
were virtually no genes (5 transcripts only) that were 
uniquely modulated by Isolate Mah 1655. In contrast, 
hundreds of genes were uniquely differentially expressed in 
response to infection with Maa 1794 and Mah VI101. 

Bioinformatics analyses demonstrated that the major 
biological functions up-regulated in response to the three 
isolates were apoptosis, inflammatory response, cytokine/ 
chemokine activity, signal transduction, cell-proliferation 
and regulation of T-cell activation (Table 2). Biological 
pathways up-regulated in the responses were related to 
apoptosis and innate immunity and cytokine signaling 
(Table 3). In addition, communication pathways between 
innate and adaptive immune responses and IL-17 signaling 
were up-regulated. Upstream regulator analysis was per- 
formed to infer the putative molecular mechanisms that 
give rise to the observed gene expression changes, and is 
illustrated in Table 4. This analysis revealed major roles for 
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Figure 3 Numbers of differentially expressed genes in human 
macrophages infected with different isolates of Mycobacterium 
avium. The Venn diagrams of numbers of up-regulated (A) and down- 
regulated (B) genes show unique and common differential expression 
of transcripts in human macrophages in response to four hours of 
incubation by three different isolates of Mycobacterium avium. Maa 
1794 is an intracellular^ replicating isolate of M. avium subsp. 
avium, whereas Mah V1 101 and Mah 1655 are isolates of M. avium 
subsp. hominissuis that only invade and persist within cells. 



TNF, NFkB, IL-1B, IFNg and others in the regulation of 
the responses. Altogether these analyses revealed a strong 
overlap between the impacts of all the three isolates. 

To identify the putative molecular mechanisms that 
are regulating the responses, we constructed molecular 
interaction networks for each isolate using mechanistic 
findings from previous studies. This analysis showed 
that Maa 1794 induced the largest gene network in hu- 
man macrophages, followed by Mah VI101, then by 
Mah 1655 (Figure 4). The networks were organized 
around a series of hyper-connected hubs, and the most 
important hubs (based on number of interactions) are 
listed in the Table 5. Notably, the network connectivity 
patterns were dominated by TNF, IL-1B, IL-6, NFkBl, 



Table 2 Analysis of biological functions induced infection 
of human monocyte-derived macrophages with M. avium 
isolates 





Biological functions 


Enrichment score 


Mah V1 1 01 compared 
to uninfected 


Negative regulation of 
apoptosis 


15.2 




Inflammatory response 


13.5 




Induction of apoptosis 


12.7 




Cytokine/chemokine activity 


8.4 




Regulation of cytokine 
production 


7.1 




Regulation of signal 
transduction 


6.6 




Regulation of T and B cell 
activation/proliferation 


6.0 




Regu ation of protein 
cytokine secretion 


5.2 




Regu ation of 
phosphorylation 


4.8 




Response to mo ecule of 
bacterial origin 


3.8 


Mah 1655 compared 
to uninfected 


nflammatory response 


15.5 




Negative regulation of 
apoptosis 


9.8 




Induction of apoptosis 


6.6 




Cytokine/chemokine activity 


5.9 




Regulation of T and B cell 

arti\/3tinn / rimlifpratirin 
dL_ 11 v a UvJ I ]/ uivJIIICIaLILJI 1 


4.5 




Regulation of 
ph os phory ation 


4.2 




Response to mo ecule of 
bacterial origin 


4.1 




Regu ation of vira 
replication 


3.6 




raiicm I ci_uy 1 1 1 uui I icLtr[JiLii 

pathway signalling 






Cell migration 


3.5 


/Waal 794 compared 
to uninfected 


nflammatory response 


16.5 




Negative regulation of 
apoptosis 


14.2 




nduction of apoptosis 


13.0 




Response to molecule of 
bacterial origin 


8.3 




Regulation of signal 
transduction 


8.2 




Cytokine/chemokine activity 


7.6 




Regulation of T and B cell 
activation/proliferation 


6.6 




Regulation of cytokine 
production 


6.0 




Regulation of protein cytokine 
secretion 


4.9 




Regulation of phosphorylation 


4.3 



The analysis was done using the Database for Annotation, Visualization and 
Integrated Discovery (DAVID). An Enrichment score > 1.3 was considered 
statistically significant. 
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Table 3 Canonical pathway analyses 





Up-regulated 


p-value 


Mah VI101 compared to uninfected 


TNFR2 signaling 


6.3E-14 




Role of macrophages, fibroblasts and endothelial cells in rheumatoid arthritis 


5.0E-11 




Death receptor signaling 


5.0E-11 




Regulation of cytokine production in macrophages and Th eel s by IL-17A and IL-17 F 


6.3E-11 




Role of IL-17A in arthritis 


1.9E-10 




TNFR1 signaling 


2.8E-10 




TREM1 signaling 


7.6E-10 




Regulation of cytokine production in intestina epithelial cells by L-17A and IL-17 F 


1 4E-09 




Hepatic fibrosis/Hepatic stellate cell activation 


1 .7E-09 




IL-6 signalinq 


1 .9E-09 


Mah 1655 compared to uninfected 


TNFR2 signaling 


6.3E-12 




4-1 BB signaling in T lymphocytes 


7.2E-10 




TNFR1 siqnalinq 


1.1E-09 




Small cell lung cancer signaling 


2.9E-09 




Role of IL-1 7A in arthritis 


4.0E-09 




Death receptor signaling 


7.6E-09 




Induction of apoptosis by HIV1 


8.9E-09 




iNOS signaling 


1 4E-08 




Type I diabetes mellitus signaling 


15E-08 




Regulation of cytokine production in macrophages and Th ce s by IL-17A and IL-17 F 


1 .9E-08 


/Waal 794 compared to uninfected 


Death receptor signaling 


1.0E-13 




TNFR2 signaling 


1.0E-12 




Role of macrophages, fibroblasts and endothelial cells in rheumatoid arthritis 


25E-11 




Small ce lung cancer signaling 


6.3E-1 1 




Role of L-1 7 F in a ergic inflammatory airway diseases 


14E-10 




TWEAK signaling 


3.0E-10 




Role of IL-1 7A in arthritis 


3.7E-10 




TNFR1 signaling 


4.1E-10 




Regulation of cytokine production in macrophages and Th ce s by IL-17A and IL-17 F 


5.1 E-10 




IL-17A Signa ing in fibrob asts 


7.1 E-10 




Down regulated 


p-value 


Mah VI101 compared to uninfected 


Role of CHK proteins in cell cycle checkpoint control 


1 .2E-04 




Tetrapyrrole biosynthesis II 


1 .9E-04 




Galactose degradation I (Leloir pathway) 


1 .9E-04 




Mismatch repair in Eukaryotes 


7.8E-04 




Heparan sulfate biosynthesis 


8.1E-04 




R/h cs f~if RRf^A 1 in P)I\IA Hsmsnp rocnrincp 
nuic ui urw^AA i n i l-'inaa uqi i lauc icjuu izat 


1 2E-03 




Heme biosynthesis II 


1.5E-03 




Glycine betaine degradation 


2.1E-03 




Coenzyme A biosynthesis 


2.2E-03 




L-serine degradation 


2.2E-03 


Mah 1655 compared to uninfected 


Colorectal cancer metastasis signaling 


4.4E-04 




Hereditary breast cancer signaling 


1.3E-03 




Role of BRCA1 in DNA damage response 


1 .9E-03 
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Table 3 Canonical pathway analyses (Continued) 



IL-12 signaling and production in macrophages 


2.2E-03 


Sulfite oxidation IV 


4.2E-03 


Antipro iterative ro e of TOB in T cell signaling 


5.1E-03 


Dendritic ce maturation 


6.0E-03 


Role of NFAT in cardiac hypertrophy 


6.5E-03 


HMGB1 signaling 


6.6E-03 


Superpathway of inosito phosphate compounds 


6.9E-03 


/Waal 794 compared to uninfected Mismatch repair in Eukaryotes 


6.3E-07 


Role of BRCA1 in DNA damage response 


2.9E-06 


CHK proteins in cell cycle checkpoint control 


5.6E-06 


CTLA4 signaling in cytotoxic T lymphocytes 


7.9E-06 


Hereditary breast cancer signaling 


1.5E-05 


Growth hormone signaling 


1.7E-05 


Molecular mechanisms of cancer 


1.7E-05 


ERK5 signaling 


2.7E-05 


Production of NO amd ROS in macrophages 


3.8E-05 


B cell receptor signaling 


1.3E-04 


Pathways induced by infection of human monocyte-derived macrophages with M. avium isolates. 
The analysis was done in Ingenuity Systems software. 
The ten top pathways are shown. 



NFkBIA, PTGS2 and others, and many of these path- 
ways were also highlighted by the upstream regulator 
analysis. 

Differential gene expression elicited by Maa 1794, Mah 
V1 101 and Mah 1655 

At FDR < 0.05, no significant differences were observed 
between cells infected with Maa 1794 and Mah VI101, 
while only two genes (IL6, GPR109A) were upregulated in 
Maa 1794 compared to Mah 1655 and one gene (IL-6) 
was upregulated in Mah VI101 compared to Mah 1655. 

Protein analysis 

Gene expression analysis pointed to apoptosis and 
immune response as two major biological functions 
upregulated in response to all three M. avium isolates, 
albeit at lower levels in Mah 1655 compared to Maa 
1794 and Mah VI101. We thus attempted to look at 
proteins involved in apoptosis by western blot and 
cytokine by a multiplex assay. However, differences be- 
tween the isolates with regards to apoptotic proteins 
were not detected, probably due to the lower sensitiv- 
ity in this assay compared to the microarray analyses 
(data not shown). 

Macrophages infected with all three isolates consist- 
ently produced more IL-10, IL-23 and TNF-oc than unin- 
fected controls (Figure 5). The production of IL-6 and 



IL-8 in uninfected cells was often high, and in some do- 
nors infection led to a decrease in cytokine production. 
The decrease in cytokine production was particularly 
evident for IL-8 after infection with Maa 1794, to a 
lesser extent by infection with Mah 1655 and rarely by 
Mah VI101 (Figure 5). There was a large donor to donor 
variation in the amount of cytokines produced, and this 
variation surpassed the difference induced by the iso- 
lates. Nevertheless, Mah VI101 consistently induced 
more IL-8, IL-23 and TNF-oc than Maa 1794 (p < 0.05) 
and more IL-23 than Mah 1655 (p < 0.05). 

Discussion 

The present study describes the uptake and intracellular 
growth of Maa and Mah isolates in human monocytes, 
and characterises genome-wide transcriptional responses 
elicited by these isolates in macrophages. Mycobacterium 
avium isolates of different genotypic characteristics were 
all able to enter and persist within human primary 
monocytes. Intracellular replication was only seen by 
Maa suggesting that the differences in the ability to sur- 
vive and replicate in macrophages cannot explain why 
humans are more often infected with Mah. The gene ex- 
pression program induced by infection with Maa and 
Mah isolates was broadly comparable in terms of the 
biological functions, pathways, and gene networks that 
were perturbed. However, the degree of perturbation 
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Table 4 Upstream regulator analyses 





Upstream 
regulator 


Activation 
z-score 


p-value of 
overlap 


Mah V1 101 compared 
to uninfected 


TNF 


12.15 


1 .40E-76 




NFkB (complex) 


9.02 


1.14E-60 




IL1B 


9.32 


9.77E-58 




CD40LG 


6.34 


2.80E-56 




TREM1 


2.56 


2.45E-49 




IFNG 


7.68 


2.73E-42 




RELA 


6.59 


3.37E-42 




TNFSF1 1 


6.16 


5.99E-41 




NFKBIA 


4.1 2 


8.41 E-40 




I KB KB 


5.4 I 


2.25E-35 


Mah 1655 compared 
to uninfected 


Nr 


8.45 


2.79E-58 




NFkB (complex) 


6.45 


1 .57E-42 




CD40LG 


4.03 


5.71 E-39 




IL1B 


6.52 


8.67E-39 




ILIA 


4.97 


1 .08E-27 




TNF5F1 1 


4.30 


1.12E-27 




RELA 


5.07 


4.74E-27 




IFNG 


5.50 


1.72E-26 




NFKBIA 


2.86 


2.75 E-26 




TREM1 


2.18 


2.80E-26 


Maa179A compared 
to uninfected 


Nr 


1 2.36 


£ one ir\ 
6.80 1-/0 




NFkB (complex) 


9.52 


6.86E-63 




GD40LG 


6.40 


2.02E-55 




IL1B 


9.70 


7.41 E-52 




IFNG 


8.72 


8.88E-44 




NFKBIA 


3.01 


1.96E-39 




RELA 


6.92 


2.09E-37 




TNFSF11 


6.55 


3.60E-37 




TLR4 


6.58 


5.10E-36 




ILIA 


7.19 


7.30E-36 



Genes predicted to act as up-stream regulators for the observed responses 
after infection of human monocyte-derived macrophages with M. avium isolates. 
The analysis was done in Ingenuity Systems software. Responses to chemicals/ 
drugs are excluded. The ten top genes are shown. 



varied considerably. The different level of expression led 
to an apparent discrepancy in the number of unique 
genes induced by the various isolates and the lack of dif- 
ferences when comparing the isolates with each other. 
Many of the genes significantly induced by M«<2l794 
were also induced at a lower level by the other isolates 
(Additional file 2: Figure SI). This intermediate expres- 
sion was neither significantly different from uninfected 



cells nor from Maa 1794. A larger samples size might 
have given more significant differences between the iso- 
lates since large donor to donor variation is present 
when using primary cell lines. 

In general, responses to Maa 1794 were excessive, and 
were characterized by increased expression of genes 
downstream of type I and II interferon signaling. Further 
studies will be required to determine if variations in 
these pathways account for the exaggerated responses 
and/or differences in prevalence. Similar pathways were 
induced in ruminant monocyte derived macrophages in- 
fected with the closely related Mycobacterium avium 
subspecies paratuberculosis (Map) [25-27], while a study 
using the murine cell line RAW264.7 showed that Maa 
and Mah induced higher levels of pro-inflammatory 
genes compared to Map [28]. Together with our data, 
these studies suggest that similar pathways were induced 
by all M. avium subspecies. A direct comparison of the 
level of macrophage activation is however challenging 
considering the variation in cell cultures and isolates 
used in the various studies. 

Macrophage activation can be interpreted as an app- 
ropriate reaction to an invading pathogen or as a failed 
attempt of the mycobacteria to silence the host's 
defences. Maa 1794 and Mah VI101 differed completely 
in their ability to replicate within host cells, although 
both isolates induced a strong activation of macro- 
phages and innate immune responses on the transcrip- 
tional level. Analyses on the protein level revealed that 
the final synthesis of inflammatory cytokines were 
higher in response to Mah VI101 compared to Maa 
1794. One could thus argue that the inflammatory re- 
sponse was sufficient to limit growth of Mah VI101, but 
not Maa 1794. The complexity of the outcome of an in- 
fection, and in vitro culture, was illustrated by the fact 
that despite a low activation of the innate immune re- 
sponse Mah 1655 did not proliferate intracellularly. The 
microarray data provided no clear explanation to the 
lower level of inflammatory cytokines induced by Maa 
1794 on the protein level compared to Mah VI101. 
However, regulation of inflammation on the post- 
transcriptional level has been demonstrated to be im- 
portant [29], and a down-regulation of immune re- 
sponses to avoid immunopathology can occur in 
chronic infections. The presence of unidentified virulence 
factors in Maa 1794 more actively shutting down the in- 
nate responses cannot be excluded. The genomes of the 
three isolates did for instance differ in the presence of 
GPL genes, and the involvement of GPLs in M. avium 
colony morphology and of ssGPLs in pathogenicity is in- 
creasingly recognised [30-32]. 

The induction of pro-apoptotic genes, such as 
RIPK2, BID and tBID was seen after infection with all 
three isolates, however apoptotic pathways were 
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Figure 4 Gene networks triggered in human macrophages infected with different isolates of Mycobacterium avium. Cells have been incubated with M. avium subsp. avium 1794 (A), A/!. 
avium subsp. hominissuis VI101 (B) and M. avium subsp. hominissuis 1655 (C) for four hours. Blue hubs indicate up-regulation, yellow indicate down-regulation of the respective gene. The size of 
the gene name reflects the number of links to other genes in the network, and genes with < 5 genes are minimized. For this analysis the number of genes was reduced, in order to comply with 
the limitation of the IPA software, by increasing the stringency of the statistical analysis to FDR < 0.01 from FDR < 0.05. 
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Table 5 Top 20 hub genes and number of connections induced by infection of human monocyte-derived macrophages 
with M. avium isolates 



Mah 1655 




Mah V101 




Maa 1794 




Gene ID 


No. of interactions 


Gene ID 


No. of interactions 


Gene ID 


No. of interactions 


TNF 


61 


TNF 


178 


TNF 


219 


IL1B 


43 


IL6 


11/ 


IL6 


141 


IL6 


38 


IL1B 


113 


IL1B 


127 


NFKBIA 


34 


NFKBIA 


81 


NFKBIA 


93 


NFKB1 


27 


IL10 


80 


FOS 


88 


ILIA 


22 


NFKB1 


70 


NFKB1 


76 


PTGS2 


19 


FOS 


69 


CSF2 


66 


NFKB2 


18 


CSF2 


59 


ILIA 


63 


CXCL2 


18 


ILIA 


58 


BCL2 


59 


CCL5 


18 


JUN 


5/ 


OSM 


59 


OSM 


1/ 


CD40 


53 


CD40 


58 


TNFAIP3 


14 


PTGS2 


52 


PTGS2 


56 


CCL3L1 


14 


Rbl 


50 


SRC 


55 


LIF 


14 


BCL2 


47 


RBI 


50 


CCL20 


14 


OSM 


47 


IL15 


50 


REL 


13 


VEGFA 


46 


CXCL2 


45 


S0D2 


13 


SRC 


43 


CCL5 


42 


NFKBIB 


11 


CCL5 


42 


PIK3R1 


42 


CCL3 


11 


CXCL2 


42 


CEBPA 


41 


TRAF1 


10 


NFKB2 


38 


TICAM1 


41 



affected to a lesser degree by Mah 1655 than by Maa 
1794 and Mah VI101. The correlation between apop- 
tosis and virulence of mycobacterial strains is debated. 
Inhibition of apoptosis has been described as a virulence 
factor in mycobacteria [33-35], but it has also been 
suggested that M. avium can use apoptosis of macro- 
phages as a tool for spreading [36]. One of the other more 
pronounced differences between the isolates was the 
increased expression of genes involved in activation of T 
cells like CD25, CD40, CD274 and IL23 in Maa 1794 and 
Mah VI101. CD25, CD40 and CD274 are transmembrane 
surface proteins involved in activation of macrophages 
and T cells upon binding, while IL-23 is essential for the 
persistence and function of Thl7 cells [37,38]. These 
findings suggest that differential activation of apoptotic 
pathways and/or activation of T cells may occur in 
humans after infection with the different isolates. The 
lower activation of macrophages by Mah 1655 may be a 
useful strategy to avoid the induction of an adaptive 
immune response in the host. This isolate was included 
because it harbours ISMpal and has a similar IS124S 
RFLP pattern to strains shown to be associated with 
virulence [24,39]. Further studies, including experimental 
infection models, could provide an answer to whether the 



"silencing approach" of Mah 1655 might be the most 
successful strategy. 

Other possible explanations for the low prevalence of 
Maa observed in the human population, and also in 
pigs, include the lack of exposure or a reduced ability to 
cross the human intestinal barrier. In the Norwegian pig 
population only Mah has been isolated [24], but the 
virulence of Maa in pigs has been demonstrated by 
others [40], and previous research demonstrated that 
Maa 1794 and Mah VI101 were equally able to infect 
piglets [41]. It thus appears that a lack of exposure is the 
most likely explanation for the difference in prevalence 
in pigs. The virulence of Maa has also been demon- 
strated in mice [42]. Although the study of Pedrosa et al. 
was performed prior to the division of M. avium into 
subspecies, one could assume that the bird isolates used 
in the older studies are equivalent with Maa, and the 
human isolates with Mah. The virulence of M. avium 
subspecies in humans can obviously not be determined 
in vivo. However, altogether the studies in human pri- 
mary cell lines and infection studies in other species, 
suggest that different exposure to Mah and Maa in 
humans is a plausible explanation for the difference in 
prevalence. 
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Figure 5 Cytokine expression from primary monocyte-derived macrophages. CD14+ cells were obtained from six donors and allowed to 
mature into macrophages by incubation for five days. The cells were then infected with the three isolates, Maa 1794, Man 1655 and Man V1 1 01 , 
at a MOI of 10:1 for 24 hours or left uninfected (neg). The supernatant was removed and assessed for IL-6, IL-8, IL-10, IL-23 and TNF-a by the 
Bio-plex assay. The six donors are represented with individual symbols. Each point is the median of triplicate wells. The horizontal line represents 
the mean response of the six donors. Mah V1 101, consistently produced more IL-8, IL-23 and TNF-a than Maa 1 794 and more IL-23 than Mah 1655. The 
difference was significant (p < 0.05) using the non-parametric Wilcoxon signed rank test. 



Conclusion 

In conclusion, our findings demonstrate that the ability 
to replicate in macrophages cannot explain the difference 
in prevalence of the two subspecies in humans. Genome- 
wide expression profiling studies further showed that simi- 
lar biological pathways and networks were elicited by the 
two subspecies, although responses to the Maa subspecies 
were exaggerated. 

Availability of supporting data 

The raw microarray data have been made publicly avail- 
able in the Array Express repository http://www.ebi.ac.uk/ 
arrayexpress/experiments/E-MTAB-1101/. The support- 
ing data are included in the additional file. 

Additional files 



Additional file 1: Table SI. List of genes differentially expressed in 
strain Maa 1794, Mah VI 101 and Mah 1655 compared to uninfected cells. 



Additional file 2: Figure SI. Volcano plots comparing gene expression 
patterns in human marophages cultured in the presence or absence of 
different isolates of Mycobacterium avium. A: Maa 1794. B: Mah VII 01. C: 
Mah 1 655. In all three plots, the red data points show the location of genes 
that were upregulated in responses to isolate Maa 1 794. 
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